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This work proposes an adaptive filter based on a new least mean sixth control 
approach with incremental conductance method of MPP for 3-phase grid- 
incorporated photovoltaic (PV) system. The proposed system comprises a PV 
array, 3-phase DC to AC converter, maximum power point tracker (MPPT), 
three-phase electronic load, and a 3-phase grid. The combination of solar PV 
array and the voltage source converter (VSC) supplies power to the grid. The 
3-phase inverter as a distribution static synchronous compensator (D- 
STATCOM) improves the quality of the system performance in case of zero 
solar irradiation. D-STATCOM also reduces total harmonic distortion (THD) 
in grid currents, improves power factor, and maintainsa constant voltage at 
the point of common coupling (PCC). The system modelling and simulation 
is achieved on MATLAB/Simulink. The proposed system performance has 
been found satisfactory and conform to IEEE-519 standards. 
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NOMENCLATURE 

Woa Wop, Wpe : Weights of fundamental active ex €p €c Adaptive active components of errors 
Components of load currents 

Wsp Wg : Total weights of fundamental active isa isp, isc Sensed grid currents 
component, Equivalent weights of 
fundamental active component 

Xpa Xpb: Xpe In phase unit templates of voltages Ka iw ise Reference grid currents 
at point of common coupling (PCC) 

Wov We : Weights of feed forward SPV Hn Adaptation constants of VSC 
power, active loss components 

Ipv» Vov Pov : PV current, PV voltage, PV power Inp Vmp»Pmp : Peak current, voltage and power of the PV 

module 

Ve : Terminal Voltage at point of FeV oe Short circuit current and open circuit 
common coupling (PCC) voltage of the PV module 

Vac» Vac “` Vaen DC-link voltage, reference value Io Rs, Rp Reverse saturation current of diode, series 


and DC-Link voltage error 


Le, R pP C F : Interfacing 


and shunt resistance of the basic PV module 
resistance 


and capacitance of RC filter 
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1. INTRODUCTION 

Due to an increase in power electronics load, three-phase electrical distribution systems are facing serious 
power quality complications such as load unbalancing, high harmonic currents, more reactive power flow, and 
poor voltage regulation. Besides, higher penetration of renewable energy further affects the power supply quality, 
and thereby, voltage levels are adversely affected. Therefore, voltage regulation has become an issue of greater 
importance in grid integration [1]-[3]. The PV array, along with the D-STATCOM converts the DC into AC, 
improves the power quality by minimizing harmonic distortions, compensates the required reactive power, and 
maintain a balance of power in all phases of electrical system [4]. It is also exchanging real power when the D- 
STATCOM is provided with an external dc source [5]-[8]. A D-STATCOM requires an efficient control technique 
for robust operation and for what various control algorithms is used, which generates pulses for switches of 
inverter. These pulses are generated through a hysteresis current controller where sensed currents are compared 
with the reference currents [9]. The principal aim of the control is to meet the load demand by taking array power 
and to supply remaining power to the 3-phase grid. If the generated solar PV power is insufficient for load, extra 
required power can be taken from the network in order to meet the load demand. In addition, control algorithms 
also offer mitigation of power quality (PQ) problems such as harmonic filtering, reactive power compensation, 
load balancing, poor voltage regulation, and power factor correction [10], [11]. 

Numerous control algorithms are used for the efficient operation of D-STATCOM, such as 
synchronous rotating frame (SRF) control, instantaneous reactive power theory (IRPT), and adaptive based 
control algorithm. The dq0 transformation based SRF control is a most common algorithm as it performs 
quite well unity power factor operation [12], [13]. However, in this control technique, the second harmonic 
component is prevailing due to load unbalances causes performance slow down. So, a low pass filter is used 
for mitigating second-order harmonic components. Other traditional controls have more response time and 
slow convergence due to the complex block, including abc-dq transformation [14], [15]. 

Adaptive filter theory exposes its power to track the atmospheric changes and characteristics of the unknown 

system. The filter parameters are self-adjusted during the changing environmental conditions that the filter 

behavior and surrounding conditions are kept to serve its purpose [16]-[19]. There are various adaptive 
controls that have been developed and available in the literature, such as recursive least square (RLS), least 
mean square (LMS). Among all of them, the LMS and its variants are highly popular for their fast 
convergence speed and stability. It is a simple and fundamental adaptive filter based control on getting the 
desired result. The main function of LMS control is to extract accurate pulse and to minimize the error for the 
efficient operation of D-STATCOM. LMS control with D-STATCOM improves current related quality 
problems such as harmonics in grid currents, load balancing of non-linear load, voltage regulation. On the 
other hand, the drawbacks of the LMS method are that when excitation is not present at the input, this 
algorithm is failed to stabilize. As a result, it produces unbounded parameter estimation and error, which put 
down the system performance and create numerical problems. Moreover, because of the low input signal, the 

LMS control algorithm also experiences stalling. Therefore, a leakage factor is added to this algorithm to 

diminish stalling while stabilizing the system at the same time and to give fast convergence response [20]-[23]. 

The proposed system, along with the least mean sixth control and single-stage topology, shown in Figure 1 (a), 

Figure 1 (b) is showing the block diagram of the proposed control algorithm showing how the reference grid 

currents are generated and then the reference grid currents and the sensed grid currents are compared using 

hysteresis controller to gererate output gate pulses for the efficient switching of the inverter. 

The proposed system provides several advantages over other topologies. 

— The least mean sixth control algorithm is a modified form of LMS algorithm. It involves higher order (fifth) 
error term in the weight updating equation. Owing to higher-order optimization, the proposed control offers 
reduced mean square error (MSE) and thereby offering reduced THD and faster convergence rate. 

— The least mean sixth control with the PV system offers lesser noise in their weights in comparison with 
the conventional LMS control. 

— The system offers greater efficiency due to one stage converter topology [13], i.e., although double stage topology 
has the advantage of flexibility in controller designing, it also has some deficiencies. With more circuit stages, the 
power loss raises, this decreases the system efficiency. Moreover, more stages increase the complexity of the 
system and thereby reduce the reliability of the system. Thus, the one-stage topology has been used. 
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The terminal voltage is calculated by equation (10) 


Calculation of adaptive error of 


P hase ‘a’, b, and jins In-phase unit templates are 
equations (12), (13), and 14) estimated by equation (11) Feed forward weight is calculate by using 
equation (17) 


Extraction of load 
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Generation of 
6 switching 
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Figure 1. These figures are; (a) general block diagram of proposed system, (b) block diagram of the proposed 
least mean sixth control algorithm 
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2. SYSTEM MODELING 
2.1. Modeling of 21 KW PV array 

Here parameters of KC200GT PV array modeled at [25] °C and 1000 W/m’ are taken as basic PV 
module [24], [25] shown in Table 1. 


Table 1. Parameters of KC200GT PV array modeled at [25] °C and 1000 Win? 


Imp 7.61 A Isc 8.61 A 
ae 26.3 V Voc 32.9 V 


o 
Pnp 200.143 W Io  9.82*10^-8 A 
Rs 415.405 ohm Rp 0.221 ohm 


Number of series units (Ns) is calculated as [13]. 


700 


N, = (=) = 7% = 26.62 = 27 module (1) 


26.3 


Number of parallel units is calculated as [13], 


N.= (Ppv/Vpv) _ (21000/700) 
p. Imp 7.61 


= 3.94 = 4 module 2) 


2.2. DC-link voltage (Vac) 
The minimum value of Vz, must be kept more than double the maximum phase voltage. So, the 
value of DC link voltage is chosen as 700 volts. 


Vac = (2V2*V,,/V3«m = (2V2 *415)/ V3 * 1 = 677.69 V 6) 


2.3. Inductor 
The value of interfacing inductance for phase leg is chosen as 2.5 mH. 


Le = (V3m = Vac)/(12jfs1;) (V3 * 1 * 700) / (12 * 1.2 * 10000 » .09 * 41.32) = 2.26 mH (4) 


Where the grid current estimated as, 


I= (= ) = —21000_ = 29.21 Amp (5) 


V3*Vn/ 1732+415 
So, peak grid current, I, = V2 * 29.21 = 41.31 Amp 


2.4. Ripple filter and DC link capacitor 

The chosen value of capacitor and resistance of the RC filter are10 uF, 5 ohm, respectively as taken 
from [13]. Here, 4600 uF as capacitance is chosen to stabilize Vdc during disturbances. The calculated value 
gives better result. 


(Pyy/Vac)/(2 * œ * Vacr) = (21000/700)/(2 * 314 *.02 + 700) = 3412.2 uF (6) 


3. CONTROL STRATEGY 

The proposed system utilizes two controls for different purposes. The incremental conductance method is 
to extract maximum irradiation from the PV system, and the adaptive based least mean sixth control is to produce 
gating pulses for the inverter. The significant equations used for the control algorithm are described as shown in. 


Ve = Tsab tsho (7) 


Vab = Asanti spe) (8) 
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Vye = Cano) 0) 
Ve = f2/3(Wsa2 + Veb? + Vee") (10) 
ee E gE (1) 
ea(n) = Ira (n)) — xpa (n) * Wpa (n) (12) 
eb (n) = Ipa(n)) — xpo (N) * Wop) (13) 
ee(n) = ILa (n)) — xpe) * Wp (n) (14) 
W.(n + 1)W.(n) + Kp{Vaen(n +1) —Veen(n) + KiVaen(n + 1) (15) 
Vaen(@) = Vac (2) — Vac(n) (16) 
Wov(n) = Zee? a7) 
Wpa(n + 1) = Wpa(n) + HnXpaCn) * (ea(n))? (18) 
Wpb (n + 1) = Wpb (n) + HnXpb (N) * (e,(n))” (19) 
Wpc(n + 1) = Woc(n) + UnXpeCn) * (e-(n)) 5 (20) 
Wg = 1/3 (Wpa + Won + Wc) (21) 
Wsp = We + We — Woy (22) 
isa” = Wep * Xpa (23) 
igo” = Wop * Xpb (24) 
isc" = Wop * Xpe (25) 


4. RESULT AND DISCUSSION 

During normal operation, the solar irradiation level has been kept at 1000 W/m? and the temperature is 
kept at 25 °C. The system remains stable during the operation shown in Figure 2 (a). Table 2 shows that the 
THD under normal operation is found satisfactory and conform to IEEE-519 standards. The solar insolation has 
been changed from 1000 W/m? to 600 W/m? at .03 second. However, the system remains stable during the 
condition. Figure 2 (b) shows the AC and DC waveforms under varying insolation conditions. Figure 2 (c) 
shows the waveforms when phase ‘a’ is disengaged through breaker for a time period from 0.03 second to 
0.06 second. During the unbalancing load period, the grid current rises, and the compensating current goes a 
little bit down, but the active power transmission from PV array to the grid remains stable. The DC link 
voltage remains constant during the operation. The waveforms of average weight component (Wg), 
fundamental weight component of phase ‘a’ (W pa), total weight component (W sp) and the feed-forward term 
(Wy) under unbalancing can be seen in Figure 2 (d). 
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(d) 
Figure 2. These figures are, (a) the system remains stable during the operation, (b) AC and DC waveforms under 


varying insolation conditions, (c) the waveforms a time period from 0.03 to 0.06 second, (d) the waveforms of 
average weight component under unbalancing 
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Table 2. The THD under normal operation 
signal Amplitude THD 
Grid voltage (Vsa) 339.6 V 1.12% 
Grid current (isa) 32.11% 3.16% 
Load current (iza) 11.25 A 28.61% 


5. COMPARITIVE STUDY 

The comparison made on the basis of the three main factors which are accuracy of extraction, signal 
oscillation and settling time. It is observed from the Figure 3 that the value of weight component (W,) 
extracted from the least mean sixth is more accurate in comparison to least mean fourth (LMF) control 
algorithm. In addition, oscillation in W, with least mean sixth is less when comparing it with the LMF 
control. Moreover, the settling time for the proposed control algorithm is less as compared to LMF control. 
The values of MSE and THD have been calculated under similar steady state condition. The MSE values for 
LMF and Least Mean Sixth cotrol are 24.9 and 17.7 respectively which means that the static error with the 
proposed new control approach is less than the LMF algorithm. Also, THD in grid currents are 4.47% and 
3.16% respectively. Table 3 displays the comparative study between LMF control and proposed adaptive 
control. 


= Table 3. The comparative study between LMF 
pom control and proposed adaptive control 
= Parameter Least Mean Proposed Control 
E 3 Fourth Control (least mean sixth) 
S Filter Adaptive type Adaptive type 
5 ; Complexity Lesser Lesser 
z ° Least Mean Sixth Degree of 4" order 6"" order 
= Least Mean Fourth rae n 
2 optimization 
g 4 Accuracy good Better 
S Static error more low 
<2 aa ops MSE 24.9 17.7 
Sampling time lys lus 
o 0.05 01 015 02 0.25 03 0.35 04 0.45 0.5 Settling time 08sec -06sec 
Time (sec) THD 4.47% 3.16% 
Oscillations less less 


Figure 3. The value of weight component (W_g) 
extracted from the least mean sixth and least mean fourth 
(LMF) control algorithm 


6. CONCLUSION 

The proposed system has been modeled for 21 KW PV power and 5.8 KW load. The system 
performance with the least mean sixth control algorithm and incremental conductance technique of MPPT 
has been used on MATLAB/Simulink. The novel control algorithm ensures rapid response to disturbances 
and fast convergence rate owing to the error increased to the fifth power. As a result, the bigger correction 
step has taken for larger errors. Moreover, the proposed system reduces the THD in grid currents and 
conforms to the IEEE-519 standard. 
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